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Abstract 

The mean field description of nuclear matter in the quark-meson coupling 
model is improved by the inclusion of exchange contributions (Fock terms). 
The inclusion of Fock terms allows us to explore the momentum dependence 
of meson-nucleon vertices and the role of pionic degrees of freedom in mat- 
ter. It is found that the Fock terms maintain the previous predictions of the 
model for the in-medium properties of the nucleon and for the nuclear incom- 
pressibility. The Fock terms significantly increase the absolute values of the 
single-particle, four-component scalar and vector potentials, a feature that is 
relevant for the spin-orbit splitting in finite nuclei. 
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1. Introduction 

The quark-meson coupling (QMC) model [0 provides a simple extension of relativistic 
many-body models based on point-like hadrons, such as quantum hadrodynamics (QHD) f2|, 
to include the explicit quark structure of the hadrons. Although both QMC and QHD and 
their variations share many common features, there exist prominent differences. Perhaps the 
most striking one is the mechanism through which the effective masses of the light vector 
mesons (p, u and 0) decrease in medium. While in QMC the decrease results from an 
increase in the lower component of the Dirac spinor of the quarks in the mesons || , in QHD 
the decrease is driven by the vacuum polarization in medium 0. Modern versions of QHD, 
based on the ideas of effective field theory incorporate vacuum polarization effects implicitly 
through phenomenological effective couplings ||. However, one can obtain either an increase 
or a decrease of the masses of the vector mesons when different sets of effective couplings 
are used - even though they fit ground-state observables of nuclei equally well. 

The QMC model has been applied to a great variety of problems in nuclear physics 
using the mean field approximation ||. In this paper we extend the model to include 
the exchange, or Fock terms, which are required by the Pauli exclusion principle. As in 
previous applications of QMC, we consider nuclear matter as a system of non-overlapping 
bags and effects due to quark-exchange between different bags are neglected. Thus the 
Pauli principle is enforced at the nucleon level. It is important to consider Fock terms not 
only as a matter of principle, but also because it is only through them that the momentum 
dependence of the meson-nucleon vertices can be explored. As we know, for example, from 
the work of Bouyssy et al. || within QHD, the momentum dependence of the exchange 
contributions from the isovector mesons (n and p) is essential to explain the magnitude and 
the systematic behaviour with N and Z of the spin-orbit splittings in finite nuclei. 

In the next section we formulate the QMC model in the language of an effective meson- 
nucleon Hamiltonian following the ideas of the Cloudy Bag Model (CBM) ||. In Section 3 
we obtain the expressions for the nuclear matter energy density and of the single-particle 
Fock potential. The numerical results are presented and discussed in Section 4. In Section 5 
we present our conclusions and discuss perspectives for future calculations. 



2 



2. Effective meson-nucleon Hamiltonian 



The Lagrangian density of the model is 



£qmc = £mit + - {d^o&o - mla 2 ) - -F^F^ + -m 2 u 
+ l - (d^ir - mlir 2 ) - -fa-GT + \m 2 p b^ 



(1) 



where £mit is the Lagrangian density of the MIT Bag Model, 





6y is one inside the bag and zero outside, 5s is a delta function on the bag surface, F^ v = 
dpUv — duUp, and G^ u = d^b^ — d v b^. The a and u;^ field operators have constant expectation 
values in symmetric nuclear matter, while 6 M field has a nonzero expectation value only in 
asymmetric nuclear matter, and the 7r field has zero expectation value in both symmetric 
and asymmetric matter, because of parity considerations. We then separate from the meson 
field operators their mean field values and treat the fluctuations in time-ordered perturbation 
theory. More specifically, we begin by solving the single-nucleon problem in the presence of 
the constant mean-fields. Then we project the quark-meson Hamiltonian, obtained from the 
fluctuating meson fields coupled to the quarks, onto the space of the single-nucleon states in 
the presence of the constant mean meson fields. 

The Hamiltonian resulting from this procedure is similar to the usual CBM Hamiltonian, 
with the difference that the effective meson-nucleon vertices in the present case are evaluated 
with density dependent wave functions. In this approach exchange effects from quarks in 
different nucleons are not taken into account. We also note that our calculation is not fully 
self-consistent, in the sense that the Fock terms are calculated perturbatively using Hartree 
self-energies. Our approach is similar to Chin's in QHD ||10|| . There, the difference between 
the perturbative and full, self-consistent results was not very large, because the Hartree terms 
provide by far the most important contributions. Since the Hartree terms also dominate in 
QMC, we expect that such a perturbative approach should be sufficiently reliable for a first 
estimate of the Fock terms in this case, and reserve for a future work a more complete, 
self-consistent calculation. 

Let o"o, c^o and &o denote the meson mean fields. Let B\(p) and B\(p) denote the creation 
and annihilation operators for single-nucleon states with spin-isospin A and momentum p in 
the presence of the mean fields. That is, the operator B\(p) creates the state 
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\p,\) = B{(p)\0), (3) 



which has energy 



E(p, A) = E*(p) + 3fl£ c + (A|r 3 |A) b , (4) 



with E*(p) = \Jp l + M* 2 . In the present paper we adopt the nonrelativistic normalization, 

(p', A'|p, A) = 5{p' - p) 5 yx . (5) 

The effective nucleon mass, M*, is obtained by solving the MIT bag equations in the presence 
of the mean fields. It is given explicitly by |JTTJ] 



= 3^L^£ + ^ B R* 3 , (6) 
R* 3 v ; 

where we have taken equal up and down quark masses, Zq is the parameter that takes into 
account zero-point and cm. motion and 



Q = yV 2 + (i?*m*) 2 , (7) 

with 

m* q = m q - g g a a . (8) 

Here, x* is the solution of the transcendental equation resulting from the linear (confining) 
boundary condition at the bag surface: 



f2 — R*m* . ^ 

Q ji(x*), (9) 



\ Q + R*m* 

where j and and ji are spherical Bessel functions. The in- medium radius, R*, is obtained 
by minimizing M* with respect to R*. 

The next step is to project the quark-meson Hamiltonian, obtained from the Lagrangian 
density given in Eq. (0), onto the single- nucleon states of Eq. (|3p. The resulting effective 
meson-nucleon Hamiltonian can be written as 

#eff = H MMF + H + W } (10) 

where -^mmf is the Hamiltonian of the meson mean fields, 

Hmmf = \™ 2 a ol - \ml u% - -mjfcg, (11) 
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H is the sum of the single-nucleon and single-meson Hamiltonians, and W is the meson- 
nucleon interaction. Specifically, Hq is the sum 

H = Y / Jdpe(p,\)Bl(p)B x (p)+ j dku 3 (k)a)(k) a] (k), (12) 

A j=cr,u),-ir,p 

where e(p, A) is given in Eq. @, ujj(k) = (k 2 + m 2 ) 1//2 , and aj-(fc) and o,j(k) are the meson 
creation and annihilation operators. The meson-nucleon interaction, W, can be written as 



W = £ / dpdp'dkS(k-p' + p) ) \ Y,Wl, x tf,p)Bltf)Bi(p)a j (k)+h.c 

j=a,LU,n,p J J2Uj[k) XX' 



(2tt) 3 / 2 

AA' 

(13) 

where the vertices Wy X (p',p) are of the general form 

W{ lx {p',p) = (p'A'|^(0)r^(0)|pA), (14) 

with 

^ = -g*, (15) 

r w = -sg (16) 

r^i^f(p'-te (17) 

r p = -^/pT. (is) 

Here = (E*(p') — E*(p),p' — p), and the polarization vectors e M (fc) satisfy 

s=l m w,p 

The next step requires evaluation of the matrix elements Wy X (p',p). For free-space 
nucleon states, matrix elements of this sort have been calculated in the MIT bag model in a 
variety of approximations. In the present case, the situation is more complicated because the 
states \p, A) have their mass-shell modified in medium. Not only is the nucleon mass modified 
from its free-space value, M, to M*, but the relation between energy and momentum is also 
modified by the mean vector fields, as can be seen in Eq. (^). These modifications are, of 
course, the interesting aspect of the QMC model as they will introduce a density dependence 
into the form factors. In order not to complicate the problem too much, in this initial study 
we adopt a simple approach. First, we parametrize the various matrix elements W{, x (p' ,p) 
as in free space, namely 
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Wl x (p',p)=u(p',\') 0*{k)u(p,\), 



(20) 



with 



0°(k) 
O w {k) 



-9lF s (k 2 



-q q e M 



0*{k)=i—G A {k 2 



-q q 



O p (k) = 

The Dirac spinors, u(p, A), are given by 
u(p, A) 



V.""(* 2 ) + 



2M* 2 



(fc 2 



T. 



(21) 
(22) 

(23) 

(24) 



J5*(p) + M* 



\l 2£*(p) 



1 

(T-p 



Xs£t 



(25) 



£?*(p) + M* 

where Xs and £ t are the Pauli spinors for spin and isospin, and E*(p) = y/p 2 + M* 2 . Next, 
we calculate the various form factors in the usual way, using, however, the bag model wave- 
functions modified by the mean fields. In calculating the form factors we ignore center-of- 
mass and recoil corrections, as well as effects due to Lorentz contraction when going to the 



Breit frame [12]. Both approximations can be improved by using the technique recently 



developed in Ref. |L3[, but no qualitative changes with respect to their density dependence 
are expected in using the present approximations. We note that in Eq. (p3|) we have not 
included the induced pseudoscalar form factor; we intend to investigate its effect on nuclear 
matter properties when using a chiral model, such as the CBM. 

Let the single-quark wave functions in the presence of the mean fields be written as 



q(r) 



g{r) 

icr-rf(r) 



(26) 



where 0(f) contains the spin-isospin wave functions. In calculating the vector meson form 
factors, we follow the standard procedure of relating the Dirac form factors F\ and F 2 to the 
Sachs form factors, Ge and Gm, as Fi = [Ge + t]Gm]/ (1 + T)) and F 2 = [Gm — + rj) 

with rj = —k 2 /AM* 2 . The various form factors in Eqs. ([H]) - fl2~4] ) are given by: 



F s (k 2 ) 
G E (k 2 ) 



drj (kr) [g 2 (r) - f 2 (r) 
dr Jo (kr)\g 2 (r)+f 2 (r) = 3G E (k 2 ), 



(27) 
(28) 
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G- M {k 2 ) = 2M* / dr 3 -^p- [2g{r)f{r)\ = G M (k 2 ), (29) 

G A {k 2 ) = drj (kr) \g\r) - ~/ 2 (r)] , (30) 

G E {k 2 ) = l -G E {k 2 ) = G E {k 2 ) 1 (31) 

G p M (k 2 ) = ~G» M (k 2 ) = ~G M (k 2 ). (32) 

Note that G E and Gm (without meson indices) are the usual, electromagnetic Sachs form 
factors. 

3. Nuclear matter energy and the Fock potential. 

The energy density of symmetric nuclear matter (in which case the mean field b does 
not contribute) is given by the sum of the mean-field energies and the Fock energy, 

£nM = 4 / 7^3 9{kF ~ E * {P) + 3 9lL °° PN + l m ° a2 ° " + ^ ^ ock ' (33) 

V / ' j = (J : UJ,7V,p 

where pn is the nucleon density and the Fock energies £p ock ,j = u,uj,Tr,p are the second- 
order exchange contributions 

Y + rrij 

Using Eqs. (pT|)-(f2"4]) in the expression for the Fock energy density, Eq. fl3j|), we can 



rewrite it as 



SU = 2 I 9{kF " A) ^(pMp, A) (35) 

with the 'Fock potential", Up ock (p), or the exchange nucleon self-energy given by 

UU = J j^e(k F -\p\)Oi{-k) ^u{p f ,X f )u(p',X')0^k), (36) 

where the four-momentum k is given as before. The potential Up ock (p) for each meson has 
the Dirac structure, 

^(P) = Ui{p) + 7° Ul{p) + T pUt. (37) 



The explicit forms of the potentials are easily obtained by making use of Eqs. (|21]) -(f24l) and 
will be given in a separate publication. 
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The final step is to determine the scalar mean field, <7 , which satisfies the self-consistency 
equation 

dp nn _ M* 1 d 

' j-a.^.K./i 



17 \ / Vr y (7 U 7=(J (J,7T,P 



where 



^ = 3^(0), C(a ) = (39) 



with 



/ s n/2 + m!i2*(n-l) , s 

sw = nin-ii'+m;^ - (40) 

The mean vector field, uo, is given in terms of the baryon density, as usual. 

4. Numerical results 

Before presenting the numerical results, we discuss some technical points. As is well 
known, the sharp surface of the bag induces oscillations in the form factors at large momenta. 
The process of removing spurious centre of mass motion and projecting onto a definite 
momentum tends to smooth this, but it is extremely time consuming. In order to avoid this 
time consuming calculation in this first investigation, we have chosen to parametrize the 
quark wave functions in terms of smooth functions - following Ref. |T|J we use a gaussian 



form with two parameters (Ro and ft). These are adjusted to reproduce the r.m.s. radius and 
the quark scalar density of the nucleon calculated with the original MIT bag wave functions. 
Note that Ro and ft are density dependent. 

We also observe that when cm. corrections and the meson cloud are neglected in the 
calculation of form factors, some physical quantities such as qa and the magnetic moments 
are not well reproduced within the model. For example, g& = Ga{0) ~ 1.09 for the "bare" 
bag without cm. corrections, whereas the experimental value is qa = 1.267 [I6|. The 
magnetic moments are given in the tensor couplings of the oj and p, 

gujfx (Uj 6 

K = lip- fin- 1 = gqpptfy = - 1 + 3 G m{0)- (41) 



Pionic corrections are known |T5[ to be very important for these quantities. For the nuclear 



matter calculation, the main effect of neglecting these is to underestimate the tensor coupling 
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of the p to the nucleoli, k v (k s is a small quantity and does not play a big role for the binding 
energy of nuclear matter). Therefore, in order to obtain a more realistic estimate of the 
contribution of the tensor coupling of the p we adjust k v = p p — p n — 1 to its experimental 
value. We also adjust GU(0) to the experimental value of g^- As we said above, the important 
aspect for our calculation is the density dependence of the form factors, and we do not expect 
qualitative changes in a more complete calculation. 

We begin by investigating the effect of adding the Fock energy to the mean field energy 
(Hartree) per nucleon. We use the a and uj coupling constants of the Hartree approximation 
to calculate the Fock terms, but solve Eq. (|38| ) self-consistently to obtain <jq and M* . The 
second term in Eq. (|T9| ) does not contribute in the case of uj coupling because of baryon 
current conservation, and for p coupling it gives an extremely small contribution to the 
nuclear binding energy and is therefore neglected. The Hartree coupling constants are fixed 
by requiring a stable minimum at E/A — M = —15.7 MeV at po — 0.15 fm -3 . For a free bag 
radius of 0.8 fm they are given |l ]| by g^/47r = 5.40 and g^/^n = 5.31, where g a is defined 
in Eq. (|39|) and g w = 3g^. We consider first the Fock terms associated with just the a and uj. 
In Figure 1 we show the results. The effect of the Fock terms is repulsive and of the order 
of 5 MeV. 

In Figure 2 we show the effect of the pion Fock term. The pion gives a repulsive con- 
tribution of the order of 8 MeV. Here, as in QHD, we have the situation that the effective 
NN interaction due to pion exchange contains a short-range contact interaction. There have 
been arguments that contact interactions should not be taken into account in a mean-field 
calculation, since they are suppressed by short-range NN correlations || . In a complete cal- 
culation, in which the full ladder diagrams are summed, this suppression of short- distance 
pion exchange is automatic. In a Hartree- Fock treatment, there is no unambiguous way to 
subtract such short range pieces from a relativistic interaction. We approach the problem 
by making a static approximation and expanding the nucleon energies in the spinors as 



The contact interaction then becomes evident and can easily be subtracted. Specifically, the 
contribution of the pion to the energy density of nuclear matter is given by 




(42) 



£ 



7T 



1 

- x 12 
2 




[G A (P'-P)] 2 
(p 1 — p) 2 + ml 




(43) 
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Expanding E* as in Eq. (f|3), we obtain 

= 6 / Wf e[kF - lpl) / Wf e(kF - |p,|) {Ga{p ' - p)|2 



2 

ml 



p' — p) 2 + m£ 



(44) 

We have recalculated the energy density using the approximate expression, Eq. ( f44|) , instead 
of Eq. (0). We found that the approximate expression induces an extremely small change 
in the result, indicating that the terms higher than 0(p 2 ) neglected in the expansion give a 
negligible contribution to the energy density. The factor 1 in the square brackets in Eq. (|44| ) 
is due to the contact interaction in the one-pion-exchange. When this factor is subtracted, 
the contribution of the pion becomes attractive (dotted line in Figure 2), and is of the order 
of 5 MeV. 

The p also gives a repulsive contribution. Because of the large value of k v , the tensor 
term gives by far the most important contribution. As in the case of the pion, the tensor 
component contains a contact term. We proceed as previously, make a static approximation 
and subtract the contact term. The energy per particle is insensitive to the static approx- 
imation. For g p (= g q p ) we use the preferred phenomenological value g p = 2 \/47r x 0.55 
and readjust the a and u coupling constants such as to refit E/A — M = —15.7 MeV at 
Po = 0.15fm -3 . The new values of the coupling constants for a quark mass of m q = 5 MeV 
and for three different values for the free bag radius, 0.6 fm, 0.8 fm and 1.0 fm are presented 
in Table I. The values of the ratios of in-medium to free nucleon masses, M*/M, and r.m.s. 
radii, r*/r, at the saturation density are also shown in Table I. In the last two columns of the 
table we present the values for the nuclear matter incompressibility and symmetry energy. 
To calculate the contribution from the Fock terms to the symmetry energy we have used the 
approximation of Bouyssy et. al, discussed on page 387 of Ref. ||. Different values for the 
quark mass between and 10 MeV do not change the qualitative results. 

Inspection of Table I reveals that the new values for the in-medium bag parameters have 
not changed much from their Hartree values ||11|| . For example, for R = 0.8 fm, we find 
that M*/M = 0.83 and r*/r = 1.02. The corresponding Hartree values are M*/M = 0.80 
and r*/r = 1.09. Also the incompressibility of nuclear matter is not changed much by the 
addition of the Fock terms. From the Hartree value, K = 280 MeV, it has changed to 
K = 285 MeV. The symmetry energy, a^, is almost entirely determined by the value of g p , 
as is well known. Our value for g p , which has not been adjusted to the symmetry energy, 
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predicts a 4 ~ 30 MeV, independently of the free bag radius. The experimental value is close 
to 33 MeV. The contribution of the Fock terms is substantial, of the order of 10 MeV. In 
Figure 3 we plot the binding energy as a function of kp corresponding to parameter set (b) 
in Table I. For comparison, we also present the Hartree solution in the same figure. 

Next we discuss the Fock potentials. It is convenient to introduce the sums V s = — g a a + 
U s and Vq = g^uio + Uq, which are the combinations that will appear in the single- 
nucleon Dirac equation. The momentum dependence of the components V s , Vq and U v at 
the saturation density are shown in Figure 4. The separate contributions of the different 
mesons to these are shown in Table II, for the three parameter sets of Table I. The addition 
of the Fock terms increases both the U s and Uq components, and the component U v is 
very small, for all three sets of parameters. For parameter set (b) of Table I, we find that 
\Vs\ + | Vq | at \p\ = kp has increased by 67 MeV, in comparison with the value in the 
Hartree approximation. Such an increase will have important consequences for the spin- 
orbit splittings in finite nuclei. However, a quantitative estimate of such effects for finite 
nuclei requires a dedicated calculation, which is complicated by the non-locality of the Fock 
interaction. We leave this for a future study. 

We have also investigated the role of the medium dependence of the form factors on 
the saturation properties of nuclear matter. As we mentioned previously, the self consistent 
change of the internal structure of the nucleon with the nuclear medium properties is the 
attractive new aspect of the QMC model, and it is important to know the role of such 
effects on our results. In order to test these effects, we re-evaluated the Fock terms using 
the form factors of free space bags, and redetermined the a and u coupling constants so as 
to obtain E/A - M = -15.7 MeV at p = 0.15 fm~ 3 . For R = 0.8 fm, the new values 
are gl/^ir = 4.37 and <y^/47r = 5.24. The ratios M*/M and r*/r remain almost the same 
as in Table I, but the incompressibility and the symmetry energy change substantially. The 
incompressibility changes from K = 285 MeV to K — 264 MeV, and the symmetry 
energy from 04 = 30 MeV to 04 = 36 MeV. These big effects arise from an increase of 
the Fock term of the a meson. For example, to the 6 MeV change in a^, the a meson alone 
contributes 5 MeV. The reason for this large effect is simple to understand. The form factor 
F s in Eq. (^) involves the difference of the squares of the upper and lower components of 
the quark wave functions, and this difference is strongly density dependent; the difference 
decreases as the density increases. This density dependence is at the heart of the saturation 
mechanism in the QMC model, and here, this effect is very visible. 
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5. Conclusions and perspectives 

In this paper we have developed a calculational scheme to introduce Fock (exchange) 
terms in the QMC model. Following the ideas of the Cloudy Bag Model, we have con- 
structed an effective meson-nucleon Hamiltonian which naturally incorporates momentum- 
and density-dependent vertices. We found that the mean-field predictions of the QMC model 
for the in-medium properties of the nucleons are maintained with the incorporation of the 
Fock terms. The nuclear matter incompressibility remains low as in the mean field approx- 
imation. Without adjusting the coupling constant of the p meson, the symmetry energy is 
predicted to be very close to its experimental value. The Fock terms increase the value of 
the four-component, single- nucleon scalar and vector potentials, which is a welcome feature 
for the phenomenology of the level splitting in finite nuclei. 

The formulation of the present model can be extended in several ways. We intend to 
incorporate the effects of chiral symmetry and the A(1232). As with nucleons in free space, 
we expect the Cloudy Bag Model to provide insight on the role of chiral symmetry on in- 
medium nucleon and nuclear matter properties. The present formulation of the QMC model 
seems to be particularly convenient for the implementation of the coupled cluster method for 
studying nucleon correlations in matter ||. We also intend to perform a fully self-consistent 
calculation of the Fock terms. With respect to the spin-orbit splitting in finite nuclei, it 
might be interesting to investigate the effect of Fock terms within an extended model with 
a density-dependent bag constant. As shown by Jennings and collaborators [T7j, when the 
bag constant is allowed to decrease as a function of the density, the description of the spin- 
orbit splittings is improved. Finally, it will be interesting to see whether, given the implicit 
density dependence of masses and coupling constants in QMC, one still needs non-linear 
meson-meson couplings for a good phenomenological description of finite nuclei. 
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TABLES 



TABLE I. a and uj coupling constants of the HARTREE+FOCK approximation, for three 
different free bag radii. Also shown are the ratios of the in-medium to free-space nucleon masses 
and r.m.s. radii at the saturation density po = 0.15 fm~ 3 , and the nuclear matter incompressibility 



K and symmetry ener 


§y a4 ■ 


The quarks mass 


is m q = 


5 MeV in all 


cases. 






R (fm) 






M*/M 


r*/r 


K (MeV) 


a 4 (MeV) 


(a) 


0.6 


4.36 


5.55 


0.82 


1.01 


289 


30 


(b) 


0.8 


4.37 


5.49 


0.83 


1.02 


285 


30 


(c) 


1.0 


4.23 


5.03 


0.84 


1.02 


273 


30 



TABLE II. Contributions of the different mesons to the single-nucleon potentials for \p\ = &F 
at the saturation density. Columns (a), (b) and (c) correspond to the parameters in Table I. The 
last line shows the HARTREE values. 



Potentials (MeV) 


V s 


V 


u v 


Parameters 




(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


Direct 




-175 


-173 


-164 


+131 


+130 


+119 








Exchange 


a 


+14 


+12 


+10 


+14 


+12 


+11 


-0 


-0 


-0 




UJ 


-44 


-39 


-31 


+23 


+20 


+16 


-2 


-2 


-2 




7T 


-4 


-3 


-3 


-3 


-3 


-3 


-3 


-2 


-2 




P 


-32 


-34 


-32 


+1 


-0 


-1 


+8 


+8 


+7 


Total 




-241 


-237 


-220 


+166 


+159 


+142 


+3 


+4 


+4 


Hartree 




-228 


-204 


-186 


+150 


+125 


+108 
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FIGURES 




FIG. 1. The energy per particle as a function of the Fermi momentum. The solid curve is the 
HARTREE result, and the dashed curve is the HARTREE+FOCK result. Here only the a and u 
mesons are included. 
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FIG. 2. The energy per particle including the a, uj and n mesons. The solid curve is the 
Hartree result, and the dashed (dotted) curve is the HARTREE+FOCK result (contact interaction 
subtracted). 
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FIG. 3. The energy per particle including the a, u, tt and p mesons. The solid curve is the 
HARTREE result, and the dashed curve is the HARTREE+FOCK result for the parameter set (b) 
on Table I. Note that the coupling constants for the HARTREE+FOCK result have been readjusted 
to produce the same saturation density and energy per nucleon. 
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FIG. 4. Single particle potentials V s = —g a <7o + U s , Vq = q^loq + Uq, and U v as a function of 
p/kF at the saturation density, po = 0.15fm -3 . The solid lines are the HARTREE+FOCK results 
for parameter set (b) in Table I, and the dashed lines are the HARTREE results. 
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